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Abstract 

We introduce a discrete time reflected scheme to solve doubly reflected Backward Stochastic Differ¬ 
ential Equations with jumps (in short DRBSDEs), driven by a Brownian motion and an independent 
compensated Poisson process. As in [5], we approximate the Brownian motion and the Poisson process 
by two random walks, but contrary to this paper, we discretize directly the DRBSDE, without using a 
penalization step. This gives us a fully implementable scheme, which only depends on one parameter of 
approximation: the number of time steps n (contrary to the scheme proposed in [5], which also depends 
on the penalization parameter). We prove the convergence of the scheme, and give some numerical 
examples. 


Key words : Double barrier reflected BSDEs, Backward stochastic differential equations with jumps, 
numerical scheme. 
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1 Introduction 


Non-linear backward stochastic differential equations (BSDEs in short) have been introduced by Par- 
doux and Peng in the Brownian framework in their seminal paper m and then extended to the case of 
jumps by Tang and Li [2Tj. BSDEs appear as a useful mathematical tool in finance (hedging problems) 
and in stochastic control. Moreover, these stochastic equations provide a probabilistic representation for the 
solution of semilinear partial differential equations. BSDEs have been extended to the reflected case by El 
Karoui et al in fjj. In their setting, one of the components of the solution is forced to stay above a given 
barrier which is a continuous adapted stochastic process. The main motivation is the pricing of American 
options especially in constrained markets. The generalization to the case of two reflecting barriers has been 
carried out by Cvitanic and Karatzas in [Tj. It is well known that doubly reflected BSDEs (DRBSDEs in 
the following) are related to Dynkin games and to the pricing of Israeli options (or Game options). The 
extension to the case of reflected BSDEs with jumps and one reflecting barrier with only inaccessible jumps 
has been established by Hamadene and Ouknine m Later on, Essaky in [5] and Hamadene and Ouknine 
in m have extended these results to a right-continuous left limited (ROLL) obstacle with predictable and 
inaccessible jumps. Results concerning existence and uniqueness of the solution for doubly reflected BSDEs 
with jumps can be found in 0,0- {SI US and 0- 
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Numerical schemes for DRBSDEs driven by the Brownian motion have been proposed by Xu in [22] (see 
also El and El) and, in the Markovian framework, by Chassagneux in [2]. In this paper, we are interested 
in numerically solving DRBSDEs driven by a Brownian motion and an independent Poisson process in the 
case of RCLL obstacles with only totally inacessible jumps. More precisely, we consider equations of the 
following form: 


f (i) Y t = Zt + Jf g(s, Y s , Z s , U s )ds + (A t - A t ) - (. K T - K t ) - jf Z s dW s - jf U s dN s , 

< (ii) Vf e [0,T], < Y t < Ct a.s., (1.1) 

l (iii) / 0 T (Pt - £t)dA t = 0 a.s. and f^(( t - Y t )dK t = 0 a.s. 

{Wt :0<t<T}isa one dimensional standard Brownian motion and {N t '■= TV t — At, 0 < t < T} is a 
compensated Poisson process. Both processes are independent and they are defined on the probability space 
(QjJtjF = {J-t}o<t<T, P). The processes A and K have the role to keep the solution between the two 
obstacles £ and £. Since we consider that the jumps of the obstacles are totally inaccessible, A and K are 
continuous processes. 


In the non-reflected case, some numerical methods have been provided: in [T] , the authors propose a 
scheme for Forward-Backward SDEs based on the dynamic programming equation and in El, the authors 
propose a fully implementable scheme based on a random binomial tree. In the reflected case, a fully imple- 
mentable numerical scheme has been recently provided by Dumitrescu and Labart in [5], Their method is 
based on the approximation of the Brownian motion and the Poisson process by two random walks and on 
the approximation of the reflected BSDE by a sequence of penalized BSDEs. 


The aim of this paper is to propose an alternative scheme to [5] to solve (1.11. The scheme proposed here 
takes the following form: 


= nr 


V? 


a j — o 


i+i 1^7 
fc” > o, 


9(tj, m 


Ojkj = o, 


,+il jy^ups + at-kj, 

e; < ?// < cr- (//.'' - s c >” = (Vj - = o. 


( 1 . 2 ) 


It generalizes the scheme proposed by [223 1° the case of jumps. Compared to the scheme proposed in [5], 
the scheme proposed here —called reflected scheme in the following —is based on the direct discretization 


of (1.11. In particular, there is no penalization step. Then, this method only depends on one parameter of 


approximation (the number of time steps n), contrary to the scheme proposed in [2] (which also depends on 
the penalization parameter). We provide here an explicit reflected scheme and an implicit reflected scheme 
and we show the convergence of both schemes. We illustrate numerically the theoretical results and show 
they coincide with the ones obtained by using the penalized scheme presented in [5], for large values of the 
penalization parameter. 

The paper is organized as follows: in Section 2 we introduce notations and assumptions. In Section 
3, we precise the discrete time framework and present the numerical schemes. In Section 4 we provide the 
convergence of the schemes. Numerical examples are given in Section 5 . 


2 Notations and assumptions 


In this Section we introduce notations and assumptions. We recall the result on existence and uniqueness 
of solution to (1.11. We also introduce some assumptions on the obstacles £ and C, specific to this paper 
(Assumption 2.5). 


Let (fi,F, P) be a probability space, and V be the predictable cr-algebra on [0,T] x Q. Let W be a one¬ 
dimensional Brownian motion and TV be a Poisson process with intensity A > 0. Let F = {Ft- 0 < t < T} be 
the natural filtration associated with W and TV. 


For each T > 0, we use the following notations: 
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L 2 (Tt) is the set of J^-measurable and square integrable random variables. 


< oo. 


• H 2 is the set of real-valued predictable processes (f> such that ||</>||jp := E Jo tfdtl 

• £J(R 2 ) is the Borelian er-algebra on R 2 . 

• S 2 is the set of real-valued RCLL adapted processes (f> such that ||<^||| 2 := E(sup 0<t<T \4>t\ 2 ) < oo. 

• A 2 is the set of real-valued non decreasing RCLL predictable processes A with = 0 and E(A|.) < oo. 
Definition 2.1 (Driver, Lipschitz driver). A function g is said to be a driver if 

• g : x [0, T\ x R 3 —» M 

(u>, t, y, z, u) i ^ g(ui, t, y, z, u) is V <S> £>(R 3 ) -measurable, 

• lls'C-,0,0,0)||oo < oo. 

A driver g is called a Lipschitz driver if moreover there exists a constant C g > 0 and a bounded, non¬ 
decreasing continuous function A with A(0) = 0 such that dF^dt-a.s. , for each {s\,y\,Z\, up), (S 2 , y 2 , z 2 ,up), 

\g(u,s 1 , y 1 ,z 1 ,u 1 ) - g(u, s 2 , 2 / 2 , z 2 , u 2 )\ < A(|s 2 - si|) + C g (\yi - y 2 \ + |^i - z 2 \ + K - u 2 \). 

Definition 2.2 (Mokobodzki’s condition). Let f, f be in S 2 . There exist two nonnegative RCLL super¬ 
martingales H and H' in S 2 such that 

Vt £ [0, T], ft <H t ~ H' < 0 a.s. 


The following Theorem states existence and uniqueness of solutions to (1.11 (see for e.g. [3] Proposition 
5.1]). 

Theorem 2.3. Suppose f and f are RCLL adapted processes in S 2 such that for all t £ [0, T], f t < ft, 
Mokobodzki’s condition holds and g is a Lipschitz driver. Then, DR.BSDE (1.11 admits a unique solution 
{Y, Z, U, a) in S 2 x H 2 x H 2 x S 2 , where a := A — K, A and K in A 2 . 

Let us now introduce an additional assumption on g, which ensures the comparison theorem for BSDEs 
with jumps (see (JOJ Theorem 4.2]). The comparison theorem plays a key role in the proof of the convergence 
of the penalized scheme (see [5]), which is useful to prove the convergence of the reflected scheme (see Section 

11 . 


Assumption 2.4. A Lipschitz driver g is said to satisfy Assumption \2.J\ if the following holds : dP (g) dt a.s. 
for each (y, z, Mi, U 2 ) £ R 4 , we have 

g(t, y, z, Ml) - g(t, y, z, u 2 ) > 0{ui - u 2 ), with - 1 < 9 < 9 0 . 

We also assume the following hypothesis on the barriers. 

Assumption 2.5. f and f are ltd processes of the following form 


ft = £0 


bids - 


aidW, 


Pl-dN s 


( 2 . 1 ) 


ft — Co + / b^ds + 


r c .dW H 


fc-dNs 


( 2 . 2 ) 


where b^, (A, < 7 ^, , p^ and pare adapted RCLL processes such that there exists r > 2 and a constant 

such that E(sup s<T \bm +E(sup 

s<T |6$n +E(sup 

S<T K n + e(su P 

s<T l°i| r ) + E(sup s<T ijsfn + 

E(sup s<T |/3|| r ) < Cj ; f. We also assume fr = fr o,.s., ft < ft for all t £ [0,T]. 
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3 Discrete time framework and numerical scheme 

3.1 Discrete time framework 

For the numerical part of the paper, we adopt the framework of m and [ 2 , presented below. 


3.1.1 Random walk approximation of (W,N) 

For n € N, we introduce S := T and the regular grid (tj)j=o,...,n with step size S (i.e. t 3 := j8) to discretize 
[0,T]. In order to approximate W, we introduce the following random walk 


W£ = 0 , 

w? = V5EUi e?, 


(3.1) 


where e",e 2 ,...,e" are independent identically distributed random variables with the following symmetric 
Bernoulli law: 

P(e? = 1) = P(e? = -1) = 

To approximate N, we introduce a second random walk 


= 0 , 

i\r n — 'S^t 4 /' 5 ] n n 

iv t — 1 li > 


where 77 ", 772 ,..., 77 )) are independent and identically distributed random variables with law 


(3.2) 


P{.Vl = K n ~ 1) = 1 - P(?7" = Kn) = Kn, 

where K n = e~ xs . We assume that both sequences e",...,e" and 77 ™, 772 , 77 " are defined on the original 
probability space (Cl, T, P). The (discrete) filtration in the probability space is F” := {J 7 " : j = 0, ..., 77 ,} with 
^0 = Ok 0} and Fj = a i e 1 > ■■■■> Vi, Vj} for j = 1,..., n. 

The following result states the convergence of ( W n ,N n ) in the Ji-Skorokhod topology. We refer to [T5j 
Section 3] for more results on the convergence in probability of ^"-martingales . 

Lemma 3.1. ( ITS . Lemma3, (III)] The couple ( W n ,N n ) converges in probability to (W,N) for the J\- 
Skorokhod topology. 

We recall that the process ( n converges in probability to £ in the Ji-Skorokhod topology if there exists a 
family M> n )n€N of one-to-one random time changes from [0, T\ to [0, T] such that sup te r 0 T i \ip n (t) — 1\ - > 0 

L ’ J n—y 00 

almost surely and sup tg [ 0 T ] — £t| —- > 0 in probability. 


3.1.2 Martingale representation 

Let 7 / J+ i denote a .F” +1 -measurable random variable. As pointed out in [T5j, we need a set of three strongly 
orthogonal martingales to represent the martingale difference m,j+i := yj+i — E(t/j+i |.7 r "). We introduce a 
third martingale increment sequence {/r" = e T )rf], j = 0, • ■ • , n}. In this context there exists a unique triplet 
( Zj,Uj,Vj ) of ^"-random variables such that 


and 


m 3+1 := y j+1 - E(y j+ 1 | J)") = V5z 0 e] +1 + u 3 rf] +1 + v 3 ^ +x , 


_ E(y j + 1 ,if +1 \Ff) _ 1 

_ _ t 

V i - E((^ +1 ) 2 |^") - «n(l-«n) 


'®(2/1+177^+11*^7) i 
E(y j+X ^ +1 \F?). 


The computation of conditional expectations is done in the following way: 


(3.3) 
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Remark 3.2. (Computing the conditional expectations) Let <3> denote a function from R 2:,+2 to R. We use 
the following formula 


, e ” +1 ■■ ,r£ + 1 )|7?) • • , e?, 1 , <,•••, rtf, - 1 ) 


+ l^$ (e?) ... , e ”, !,<,••• ,m", K „) 


2 

1 — k. 


, e", — 1 , 77 ", • • • ,rjf,K n ). 


3.2 Reflected schemes 


The barriers £ and £ given in Assumption 2.5 are approximated in the following way: for all A; E {1, • • • , n} 

fc-i 


k -1 


/c — 1 


4 n = 4 + E+E 4^r +1 + E &*>?+ 1- 

4=0 4=0 4=0 

Ac—1 fc-1 /c-1 

cr = Co + E^+E 4^ e ? +1 + E 4<+i- 


(3.4) 


(3.5) 


4=0 


4=0 


4=0 


Lemma 3.3. Under Assumption 2.5 , there exists a constant depending on C^ t T and X such that 

^supmaxE(|^| r ) + supmaxE(|C”r) + supE(|^| r ) +supE(|^| r ) < C s f;T;A 

n J n 3 t<T t<T 

(ii) £" (resp. ( n ) converges in probability to £ (resp. () in Ji-Skorokhod topology. 

Proof, (i) ensues from Burkholder-Davis-Gundy and Rosenthal inequalities, and (ii) ensues from jHi, The¬ 
orem 6.22 and Corollary 6.29]. □ 

In the following Section we introduce the implicit reflected scheme, which is an intermediate scheme 


useful to prove the convergence of the reflected scheme ( 1 . 21 . 


3.2.1 Implicit reflected scheme 

After the discretization of the time interval, our discrete reflected BSDEs with two RCLL barriers on small 
interval [tj,tj+ 1 [, for 0 < j < n — 1 is 


E = Vj+1 + Zf, ups + a n - k? - z?V 6 e ? +1 - uffl +1 - vfyf +1 , 

af > 0, kf > 0, afkf = 0 , G <yf < C, (yf - )a? = (yf - G)kf = 0. 


(3.6) 


with terminal condition y™ = fff. By taking the conditional expectation in (3.6 1 w.r.t. J-f, we get 


' y n = C 

on z>n 


(s ) J v? = %?+il+ g(tj,yp q,ups + a™ - k? 

' 1' l n. m \ r\ n. l n r\ 


a] > 0 , k™ > 0 , oijk'j = 0 , 

,c; < //; < c;. (yj - qw; = (y? - cpis; = o. 

Lemma 3.4. For S small enough, (5i) is equivalent to 

r vi = c, 

(S J »" = + a i - k 3 ). 

2J | a" = (Elyp^p + g(tj,q, zp upS - q)~, 
U" = (nyp i \Tp+g(t J ,q,zpupS-Cp + , 

where T(y) := y - g(tj,y, Zj,upS. 
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Proof. For 5 small enough, T is invertible because the Lipschitz property of g leads to — ^(y'))(y—y') > 
(1 - 6C g )(y - y') 2 > 0 for any y ± y'. 

We first prove that (Si) implies (S 2 ). Let us firstly assume that Vj <n — 1, £" < On the set {yj = £"} 
we have k? = 0, then af = 'F(^) - E[^ +1 |J"«] = (E[y? +1 \F?} - tf(Q*))- (since E[y? +1 \jg] - *(#) = 
^(y") - *(f?) -af < 0) and on {yf > £"} we have af = 0, (E[y? +1 | F?] - T(^) = ’%?) - * (£?) + kf > 0 
(thanks to the monotonicity of 4')). Then, a" = — \l/(£")) _ . The same type of proof leads to the 

fourth line of (S 2 ). If there exists j < n— 1 such that £" = (f, we get £" = £" = yf. Then, we have o" = 0 or 
kf = 0. If both are null, we get \F(t/") = E[y” +1 |J r ”] = \F(£”) = 4'(£J ! '). This coincides with the definitions of 
af and kf given in (S 2 ). If a? > 0, kf = 0 and we get af = - E[y? +1 \J?] = 4 f (^”) — E[j/” +1 1 J"”], then 

a" = (E[y? +1 |J7]-*(£?))-. Conversely, assume (S 2 ), let us prove afkf = 0, iy'f - )a r f = (//” q)k’‘ = 0 

and £" < yf < If a" > 0, we get \F(£") > \F(£") > E[yf +1 \Fj], then kf = 0. Let us prove that 
(Vj — £”)a" = 0. If a’j > 0, 4'(i/") = E[j/” +1 |J 7 ”] + a" = \F(£"). Since T is a one to one map, we get 
yf = The same argument holds to prove (yf — C, ”)fc" = 0. Let us prove that £" < yf. To do so, 
assume that yf < In this case a" = fc” = 0, which gives \F(£") < E[y" +1 |JJ l ], by definition of a". Then 
T (yf) = E[j/" +1 |J^] > \F(£"). 4> being a non decreasing function, this leads to absurdity. □ 


We also introduce the continuous time version (Y) n , Zf, Uf 


Wr,AT)o<t<T Of (jtf 


i l,n\ 
i 1 K j ) 


j<n : 


[t/8] [t/8] 

V n — n n 7™ — v n 77™ — ,,™ 4™ — \ ' n n jyn — \ rn 

*t ■—y[t/s\i Zj t ■— z [t/8]i u t -— u [t/8}i A t •— ~~ / ; • 

2—0 2—0 


(3.7) 


In the following 0 n := (F n , Z n ,U n ,A n - K n ). 


3.2.2 Explicit reflected scheme 

The explicit reflected scheme is introduced by replacing yff by E[?;” +1 |J r ”] in g. We obtain 


[y” = y] +} +g(t J ,E[y^ +1 \Tf ] 1 z],u-)8 + a] - k] -u ^ +1 

< a] > 0 , fc" > 0 , a?fc” = 0 , 

U" < v? < g, (57 - gpf = (yf - qji% - 0 . 


L IL 

j A 6 jF h- 15 


(3.8) 


with terminal condition y" 
derive that: 


£". By taking the conditional expectation in (3.8 1 with respect to J 7 ", we 


(5i) < »7 = E[S7 +1 |J7] + 5 fe,E[y7 +1 |J 7 ;],^ ! u|)(5+ 07 - fc” 57 > 0, k ] > 0, = 0, 


'Vn=tLn, 

y? = m 

[g < //; < c”, (7/; - gp? = (v? - Q) k l = "• 

As for the implicit reflected scheme, we get that (5i) is equivalent to (6> 2 ) 

r y n = p 


We 


<£S| 

^■3 

II 

nv?+i\??] 

+ 

E[y? +1 |J7], 

0",tZ7)(5 + 


. T” 

> 

= 

(E[^ +1 |J7; 

\ + 9(tj 

>E[y"+i -A"] 

,zf,u?)6- 


5 

?$-| 
^S. S 

II 

(E[y7 +1 |J7 

} + g(tj 

.E[y7 +1 |^7] 

,z?,*Q)5- 

-CD 


- TL - TL - Tl - TL - Ti 

continuous time version (Y f ,Z t ,U t ,A t , K t )o<t<T of (yf, z‘ 





[t/8] 


[t/5] 

= V[t/8 ]> 

y n - ~n 

■— z [t/8\ 

u n — 

5 u t •— 

u [t/S]’ A t •- 

V<-A 

2 — 0 

n 

t : = 

E^- 

2=0 


7 —77 7 

i, a j , kj )j< n : 


(3.9) 


In the following 0™ := (h", Z”, {/”, A" — /\") and a” := A" — Ai". 
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3.3 Implicit penalization scheme 

In this Section we recall the implicit penalization scheme introduced in [5]. The penalization is represented 
by the parameter p. As the implicit reflected scheme, this scheme will be useful to prove the convergence of 
the explicit reflected scheme. For all j in {0, • • • , n — 1} we have 


'yf n = 1 + zf n , uf n )6 + a p f n ~ kj’ n - (z? n V6e? +1 + uf n V ? +1 + v? n f? +1 ), 

= P 5(yf n - $)~, fcf" = p6(Q yf n )~, 

v P,n cn 
. on ’ Sn' 


(3.10) 


Following (3.3), the triplet (z?’ n ,iiE’ n ,v I l' n ) can be computed as follows 


^ = 1*7). 


Taking the conditional expectation w.r.t. JFf in (|3.10 1 , we get 


2 /r=( ^ ) -i (E(2/ ^ i|7?))) 


3 


= vKvT - c;r- K' n = p s W i/T)“ > 


1 


where = y - g(jS,y,zf n ,uf n )6 - p6(y - + pS(Q-y) . 

We also introduce the continuous time version (Yf' n , Zf’ , Uf ! ' n , A p,n , /\f’”)o<t<T of the solution of the 


discrete equation (3.10|): 


[*/<*] [t/S\ 

vP,n _ p,n yp,n _ p,n Tr p,n _ p,n >p,n _ p,n t/P,™ _ U>,n 

*t — V[t/5}’*t ■— Z [t/5\' U t ■~ U [t/S]^ /i t ■— a i ’ 

i—0 i—0 

and a p,n := A p ’ n — K p ’ n . The following result ensues from J5j Theorem 4.1 and Proposition 4.2]. 


(3.11) 


Theorem 3.5. Assume that Assumption \%.5\ ho lds and g is a Lipschitz driver satisfying Assumptio n \ 2.f\ 
The sequence (Y p ' n , Z p,n , U p,n ) defined by (3.111 converges to (Y,Z,U), the solution of the DRBSDE ( |1.1[ )' 
in the following sense: Vr £ [1,2[ 

lim lim ( E 

p —»oo n—t oo \ 

Moreover, Z p,n (resp. U p,n ) weakly converges in H 2 to Z (resp. to U) and for 0 < t < T, converges 

weakly to at in L 2 (Ft) as n —> oo and p —> oo, where ( if n ) n eN is a one-to-one random map from [0,T] to 
[0,T] such that sup tg r 0 T i \if n {t) — t\ - > 0 a.s.. 

L ’ J ti .— 


[ \Yt n ~Y s \ 2 ds 

+ E 

[")Z*> n -Z s \ r ds 

+ E 

f T \ur-UsU3 

) = 0 . 

(3.12) 

Jo 


Jo 


Jo 

J 



4 Convergence result 

We prove in this Section that 0" converges to 0 := (Y t , Z t , U t , A t — K t ) 0 < t < T , the solution to the DRBSDE 


(1.1). The main result is stated in the following Theorem. 


Theorem 4.1. Suppose that Assumption \2.5\ holds and g is a Lipschitz driver satisfying Assumption \2.f\ 
Then we have 


lim E 

n—¥oo 


[ \Yt-Y t \ 2 dt+[ \Z n t -Z t \ 2 dt+ [ \U n t -U t \ 2 dt 
I o Jo Jo 


= 0 . 
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Moreover, oQn/ t \ converges weakly to a* in L 2 (Ft). 


Proof. To prove this result, we split the error in three terms. The first one is the error 0 — @ n , the second 
one is 0™ — 0 p,n , where 0 p,n ■= (Y p ’ n , Z p ’ n , U p,n , A p,n — K p ’ n ) represents the solution given by the implicit 
penalization scheme (see (3.111), and the third error term is 0 p,rl — 0, whose convergence has already been 


proved in [5]. The result on the convergence of 0 p,ra to 0 is recalled in Theorem 3.5 


U) 


We have the following inequality for the error on Y (the same inequality holds for the errors on Z and 


E[ / | Y t - Y t \ 2 dt] < 3E[ / | Y t - Y t n \ 2 dt\ + 3E[ / \Y t n - Y t p ’ n \ 2 dt} + 3[ / \Y p ’ n - Y t \ 2 dt\. 


For the increasing processes, we have: 

mut) ^ a *i 2 ] ^ 3 - ^ (t) i 2 ]+EK» (t) - 

Then, combining Propositions |4.5] |4.6|and Theorem |3 . 5| yields the result. 


r\ 2 ]+n\ar-a t f. 


(4.1) 

□ 


Definition 4.2 (Definition of c and No). In this Section and in the Appendix, c denotes a generic constant 
depending on C g , ||< 7 (-,0,0,0)Hoc andC^^ t \,T- N 0 is defined by N 0 := 4T(1 + C g + C 2 + ). 

The rest of the Section is organized as follows: Section |3.3| recalls the implicit penalization scheme 
introduced in [5] and the convergence of 0 p,n — 0, we g ive some intermediate results in Section 4.1 and we 
prove the convergence of 0 — 0” (see Proposition 4.51 and the convergence of 0" — Q p n (see Proposition 


4.6) in Section 4.2 


4.1 Intermediate results 

In this Section we state two intermediate results useful for Section 14.21 
Lemma 4.3. Under Assumption WT^we have 


supE[|y: 




-E 


*E 

i=o 


Z?\ 2 + K; 


lb — X lb — X 

.(1 - «„) Y. Kf + V K 

3=0 j =0 


n— 1 


5 E K\ 2 

3=0 


< c. 


Proof. Since < yf < C\ Assumption 2.5 gives sup ■ E(|y”| 2 ) < c. Let us deal with zf and u r j. To do this, 


we apply Lemma |B.1| with io = i and i\ = * + 1 to the process y n and we sum the equality from i — j to 
i = n. We get: 


E II^| 2 ]+^ E £ [l-n 2 ] + «n(l - «n) E E tl<| 2 ] 


l=J 


l=J 


<E[|C| 2 ] + 2(5 Y, nVi9{tu V?, Z?, <)] + 2 £ VO ~ 2 E VO. 


i=j 

n—l 


i=j 


i=j 


— E [|?n | 2 ] + ^ E 0,0,0) 2 + S I 1 + 2 C g + 2C 2 + 


2C 2 5 


t=j 


^n(l ^n) 


E E H y r i 

i=j 
n—1 


n |2i 


f ^ 1 /-i \ 7T. 1 q r 

+ \ E E [k"i 2 ] + M 9 Kn E E [i<i 2 ] + -E E (iy"i 2 ) + ^E E (Ki 2 ) + j E E (i fe ?i 2 )- 


i=3 


%=3 


i=3 


i=3 


%=3 


2 













Since < y™ < C", we get 

< < (E(tf +1 |S?) + ^(hkkMO - £*)" = Kbi +9(ti, 

k? < (E(cr+1 \Q?) + Sg(M, C, M <) - cr) + = S(bi + g(U, Q, z?, <))+. 

Then, using the Lipschitz property of g gives 

n —1 n —1 

f ^E(|a”| 2 )<5a<5^E[|6«| 2 +|^,0,0,0)| 2 + C' g 2 (|C| 2 +k"| 2 +l<r)] ) 


(4.2) 


(4.3) 


and the same result holds for f J2i=j E (|fc?| 2 )- By Using Assumption 2.5 and the inequality sup, ; E(|y"| 2 ) < 
c, we get 


n— 1 


n— 1 


« E E [N"| 2 ] + Ml - *») E E tl<| 2 ] < C + 5 M 10aC 9 ) E E (N"| 2 ) 


n— 1 


+ k „(1 — n n ) ( — + lOaC” 


K n (l - K n ) 


n— 1 


E E (i<i 2 )- 


Since (1 _k) Ktl = j (1 _ e -k) e ~ A<i and e x < fkq < e 2x , we get (1 _ K d „) Krl < {e 2X1 . Then, by taking a = 
40^2 (Ae _2AT A 1), we get 5 YJiZj E [k"| 2 ] + M 1 - K n ) YJiZj E [|<| 2 ] < c. Plugging this result in (JO) ends 
the proof. 


The same type of proof gives the following Lemma 
Lemma 4.4. Under Assumption\2. -7[ we have 


su P E[|y"f]+E 


5 E 

3 =0 


^r+«, 


.(i-gEi“"i 2 + 


E 


S?l 2 


j-0 j =0 

4.2 Proof of the convergence of @" — @ n and 0" — 0 p,i 


1 n — 1 

.Eik n i 


J=0 


< c. 


Proposition 4.5. Assume that Assumption 2.5 holds and g is a Lipschitz driver. We have 


□ 


lim sup E[|r; - 171 2 ] + E[ / \Z‘ - Z:\ 2 ds\ + E[ / \U S - U:\ 2 ds} = 0. 

n-*-°°o< t <T Jo JO 


(4.4) 


Moreover, limn-^(a" — a") = 0 in L 2 (tF t ), for t £ [0,T]. 


Proof. Let us consider yf , the solution of the discrete implicit reflected sheme (3.61 and yf, the solution of 
the explicit reflected scheme ( |3.8| ). We compute yf — y 7 j | 2 , we take the expectation and we get: 

n\y] v] I 2 ] <E[IMi - Mil 2 ] - ®0k" - k"l 2 ] - Mi - M E 1K - M 2 ] 

+ 2«5E[(^-y70(t j ,y77, U 7- 5 (t J ,E[y7 1 |^],77u7)] 


-E 


«(s(ij t »?>Mk) -sk, E lMiiMMkk")) + K -k") - (k - k”) 


n 2 


+ 2E[(y" - y])(a] d])} 2E[(y? - y])(k? - fc”)], 

< E[|j/” +1 - ^ +1 | 2 ] - - 3?| 2 ] - K „(l - «„)E[|«7 - k>| 2 

+ 25E[(y" - yk0k,2/k^kk") ~ 3k, E [MM"MkM>)]' 
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The last inequality comes from (y™ — y")(a" — a”) < 0 and (yf - — /cj) > 0 (this ensues from the 

third and fourth lines of ( 1 S 1 ) and (5i)). Taking the sum from j = i to n — 1 we get 


n— 1 


n—1 


E[k" - v?\ 2 ]+s £ e[|^ - *£] + Kn (i - «„) £ e[|u” - u£] 


n— 1 

< 2<5 £ E[(y? - yj)(g(t j ,y'j, z?, v?) - g{t jt E[y^ +1 \F?],z?,u?))], 

n —1 ✓ X \ n— 1 

< 2 SO, £ E [I#? - V?IK - E[5? +1 l^]l] + 2 SCI (l + £ E[l»? - 5?t] 

J- 4 J - 1 


+ 5 E e 1|2"-' 


i (1 ^n) 


£e[|«?- 


^ n l 2 ] 


(4.5) 




J=l 


Since y” - E[y? +1 1JJ 1 ] = y] - y] + y] - E[y? +1 \J*] = y? - y? + Sg{tj, E[j£ +1 IJ?], 2?, u?) + S? - *£, we get 

2<5C S E [|y? - y"||y" - E[^ +1 k;]|] < (2C S + l)ffi[|^ - y£] + C 2 ffi [(|5 5 (t J -,E[|^ +1 |75*],^,«y)| + Kl + 


Plugging the previous inequality in (4.51 and using Lemma 4.4 gives 

n—1 


r n — 1 /i \ n —1 

E[|y" - V? I 2 ] + 2 £E[|^ - 2?| 2 ]+ M ~ Kn) £ E[| 

j=i 

( . _ ... 2C 2 5 




n—1 


< I 1 + 2C S + 2C 2 + 

\ 9 Kn(l-Kn) 


<5£E[|y”-^| 2 ] + c5 2 . 


Let n be bigger than N 0l then <5 ^1 + 2 C g + 2C 2 + K ^ ^ < 1 (for all n > 1 we have K ^ < ^-e 2AT ). 
The assumption on S enables to apply Gronwall’s Lemma to get sup 0<i< „E[|t/" — y”| 2 ] < cS 2 . Plugging 


this result in the previous inequality leads to (4.4 1 . The convergence of ( A n — K n ) — (A™ — K") ensues from 


A n t - K? = F 0 ” - Y t n - 


K - K r ; = y£ y" - 


g(s, y s n , z:, u:)ds + T zyw: + f u?<w?, 

Jo Jo 


I g(s,Y n s ,Z n s ,U n s )ds+ I Z,dW? + I U'fdN ]\ 
'o Jo Jo 


from the Lipschitz property of g and from (4.41. 


□ 


Proposition 4.6. Assume that Assumption 2.5 holds and g is a Lipschitz driver. For n > No, we get 


sup E[| Y t n - Y t p ’ n \ 2 ] + E[ [ T \z: - ZP’ n \ 2 ds\ + E[ [ T I u: - ur\ 2 ds\ < £. 

0 <t<T Jo Jo \fp 

Moreover, V t £ [0,T], E[|a" — af ,Tl | 2 ] < ■ 


(4.6) 


Proof. Let us first prove (4.6 1 . From (3.61, (3.10) and Lemma B.l applied to the process (y n — y p,n ) following 
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the beginning of the proof of Lemma |4.3| we get 

n— 1 n— 1 n— 1 

m - yf n i 2 +^ E E i- n - *n 2 + a - «»)«„ e nw? «ri 2 ] + a - «„k e E ik n - <’"i 2 ] 


»=7 




*=7 




*=J 
n—1 


+2 e e[(w? - yn w - of’")] - 2 e E i(yr - »f' n )(*r - *r n )]- 




4=7 


Let us deal with the last two terms 


(y? - vDW - <’") = (v? - &M - (yf’" - CW - (y? - $K n + (yf" - £>f" < (yf" C)“ <■ 

By using same computations, we derive 

(: y r - yf")(fc? - *f") > -(yf" - (?) + K- 

By using the Lipschitz property of g , we get 

n\y^-y? n \ 2 ]+\mq *ri 2 i+-«ri 2 i 


/ 2 2Cg(5 

< 2C S + 2Cg + 9 

\ ^n(l l^n) 


n— 1 


n—1 


<5 E E[(y? - yf ") 2 ] + 2 E E I(yf n - C)"o? + (yf" - CD**?]- 


Using Cauchy-Schwarz inequality gives 


E[|y" - yf"] 2 ] + Le[|z? - zf ”| 2 ] + Wn) E[|^ - uf " | 2 ] 


< 


2C?<S 


2C ” +2C » + K„(iU„) J 15 E E i« - »•’") 2 i 


2 C*5 


i=3 



1 7t— i \ / 71 — i 1 7t—i 

- E £[(«?) 2 ] + 2 [8 E E [((yf’ n - C) + ) 2 ] t E E K*") 2 ] 


*=? 


4 =< 7 





l=J 


71 — 1 


< 1 2c 9 +2C 9 + „ ;r g f,, ) <5 e e m - yf") 2 ] 


^n(l ^n) 


n—1 


+ ^(sg*K«r)i) ig^iw) 2 ] 


i—j 

n—1 

(5 


^(^EW”) 2 ! ^E E [«) 2 ] 


i=J 



iV. 


»=7 




Since n > iVo, Lemma 4.3 Lemma A.l and Gronwall inequality give (4.6). Concerning a™ — a^' n we have 

a? - «r =(1? - 5?’") - (F 0 n - Y 0 p ' n ) - f g(s, Y s n , Z?, U s n ) - g(s, Yf’ n , Z?’ n , U^ n )ds 

Jo 

+ f (z? - zrww: + [ (u: - ur)dN ?. 

Jo Jo 

It remains to take the square of both sides, then the expectation, and to use the Lipschitz property of g 
combining with (4.6 1 to get the result. 

□ 
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5 Numerical simulations 


We consider the simulation of the solution of a DRBSDE with obstacles and driver of the following form: 
:= {Wt) +2(1— ^)N t + |(T— t), (t '■= (Wt) + (1— yO((N) 2 + l) + \{T — t), g(t, u, y, z, u) := — 5\y+z\+6u. 

Table [I] gives the values of Y 0 with respect to n. We notice that the algorithm converges quite fast in n. 
Moreover, the computational time is low. 


Table 1: The solution y n at time t = 0 


n 

10 

20 

50 

100 

200 

300 

400 

Vo 

1.2191 

1.3238 

1.3953 

1.4167 

1.4293 

1.4332 

1.4352 

CPU time 

2.14 x 10" 4 

1.5 x 10" a 

0.0211 

0.1622 

1.4230 

5.2770 

12.5635 


When we use the explicit penalized scheme introduced in 0> we get yf/' = 1.4353 for n = 400 and 
p = 20000. The CPU time is 12.85s. 

Figures [l] [2] and [3] represent one path the Brownian motion, one path of the compensated Poisson process 
(with A = 5) and the corresponding path of (yf, Cf)i<i<n- We notice that for all i, y™ stays between the 
two obstacles. The values of and y^’ n are almost the same when n = 400 and p = 20000. The CPU times 
are also of the same order. The main advantage of the reflected scheme is that there is only one parameter 
to tune (n). 


one Brownian path 



Figure 1: One path of the Brownian motion for n = 400. 


A Technical result for the implicit penalized scheme 

In this Section, we use N 0 and c introduced in Definition |4.2| 

Lemma A.l. Suppose Assumption \2.5\ holds and g is a Lipschitz driver. For each p £ N and n > N 0 we 
have 

n— 1 n —1 1 n —1 1 n —1 

supE[|#"| 2 ] + 5 £ E[| 2 ff] + K n (l Kn ) J2 E[K’"| 2 ] + ^E E[i«ri 2 ] + i E E[|*TI 2 ] < c - 

4 j —0 j —0 j —o j=o 

Proof. By applying Lemma |B.1| to the process y p ’ n between i and i + 1 and by suming the equality from 
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compensated Poisson path 



Figure 2: One path of the compensated Poisson process for A = 5 and n = 400. 


trajectories of the solutions and the barriers 



Figure 3: Trajectories of the solution y n and the barriers and for A = 5 and n = 400. 
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i = i to i = n, we get 


n— 1 


n— 1 


n—1 


E[|#T] + 5 E E D*fT] + «n(! - «™) E E [K’"| 2 ] + Mi - «n) E 

i=j i=j i=j 

n —1 

< E[|C! 2 ] + + 2EE(sTX" " tf’"*?’")]- 


n—1 


Note that yf’X’" = -^K’") 2 + " and = W n ) 2 + «’"■ We have that: 


£[| 2 / : 


P,wi2i 


r n—1 

D 


p,ni2i 


n —1 


i (1 ) 


n—1 

E> 


p,n |2i 


1 n— 1 

^E ] 


n—1 


a: 


P,ni2i 




p<5 


< E[ic:| 2 ] + <5£[E 1 9(U,0, 0,0) I 2 ] + 2<5 1 + 2C g + 2 C 2 g + 


2CiS 


1=0 

2 c- \ n—1 


Kn(l - «n) 


E E ^r 


n—1 


n—1 


+ 2 £ E[(er)aH - 2 E EKCDfcH- 


We get 2EE E [(a«H < aE(su Pi |ff| 2 ) + ^E (e"” 1 a?’")* and 2 E?7 E[(0*)fc?'"] < /?E(su Pi |Cf) 
dfi (E'E • Following the same type of proof as [TB] Lemma 2], we get 

/ n —1 \ /n—1 \ n —1 

e + E E*? ,n <c(c+EE«(ri^i^iv^^)wi^i^ 






*=? 


Finally, by taking a = (3 = 4(7 and by applying the Gronwall inequality (we recall n > Nq), we get that: 


supE[|^’ n | 2 + T > \z: 


lE 

J=0 


,n|2 i ^n(l ^n) \ ' | p,n|2 , 1 \ ' 

+ 4 I + La 


a 


p,n |2 


J=0 


p<5 ' 1 0 

y j=0 


s E i*?ti < «=• 


pS z 
* o=o 


□ 


B Some results on discrete stochastic calculus 

In this section we present two lemmas which are used throughout the paper. 

Lemma B.l. Consider two integers io and i\ in 0, ...,1V and (y n ) n a discrete process. We have 

Vi i = vl o + 2 VioiVn - Vi 0 ) + (Vn - Vi 0 ) 2 ■ 

The proof comes from the computation of ((& — a) + a) 2 , we omit it. 

Lemma B.2. (A discrete Gronwall lemma) Let a, b and a be positive constants, Sb < 1 and a sequence 
( Vj)j = i t .._ n of positive numbers such that for every j 

o 

Vj + a < a + bS Vi. 

i=l 

Then 

bT 

sup Vj + a < ae . 

j<n 

A proof of this lemma can be found in m Lemma 2.2], so we omit it. 
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